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The aging of operational reactors leads increased mechanical vibrations of reactor internals. The vibration 
of the in-core sensors near their nominal locations is a new issue for the neutronic fields reconstruction. Cur- 
rent field reconstruction methods fail to handle spatially moving sensors. In this work, we proposed a Voronoi 
tessellation techinque in combination with convolutional neural networks (V-CNN) to handle this challenge. 
The observations from movable in-core sensors are projected to the same global field structure, this projection 
is achieved with Voronoi tessellation, holding the magnitude and location information of sensors. The general 
convolutional neural networks were used to learn the map from observations to the global field. The proposed 
method is able to reconstruct the multi-physics fields (e.g., the fast flux, thermal flux and power rate) using 
observations from single field (e.g., thermal flux). Numerical tests based on IAEA benchmark proved its po- 
tential for real engineering usage, particularly, within an amplitude of 5 cm around nominal locations, the field 
reconstruction leads to average relative errors below 5% and 10% in L2 norm and Lo norm, respectively. 
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I. INTRODUCTION 


Since its advent in the 1950s, nuclear energy plays a cru- 
cial role in meeting world energy needs and is also an im- 
portant component of clean energy today. Nuclear energy is 
mainly generated through nuclear reactors, which are gener- 
ally designed to operate for 30 to 40 years and can last even 
longer with license renewals. Based on data from the Power 
Reactor Information System (PRIS), among the total of 437 
reactors, 289 reactors have been in operation for more than 
30 years [1]. In other words, more than 6046 of the current 
nuclear reactors are facing aging issues. As a consequence, 
reactor operational problems or anomalies are expected to be 
more frequent. The aging of operational reactors also leads to 
increased mechanical vibrations of reactor internals such as 
core barrel, control rods, in-core instruments and more specif- 
ically fuel assemblies, or other vibrations such as flow block- 
age, coolant inlet perturbations [2—6]. 

Various reactor core monitoring techniques aim to address 
these challenges, mainly based on using the observations of 
neutron flux acquired by in-core and ex-core instrumentation, 
combining numerical simulations. Such techniques and sys- 
tems include, but are not limited to CORTEX [7], BEACON 
[8], RAINBOW [9], etc. We refer the readers to [10] a detail 
overview of the reactor core monitoring techniques. The pro- 
cess that combining observed data in the core and the simula- 
tion data is generally called field reconstruction [11, 12], with 
the main goal to infer the neutronic field in the core, thereafter 
the safety related parameters can be calculated, such as en- 
thalpy rise hot channel factor (FdH), peak heat flux hot chan- 
nel factor (FQ), linear power density of fuel rods (LPD) and 
deviation from nucleate boiling ratio (DNBR), etc. 

On the key algorithms for field reconstruction, one proto- 
type is data assimilation that originally arises in earth sciences 
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including meteorology and oceanography [13]. Data assim- 
ilation framework allows combining observations and mod- 
els in an optimal and consistent way, including information 
about their uncertainties [14—16]. It has been applied in sev- 
eral studies in nuclear engineering fields [9, 17—21] for field 
reconstruction in a unified formalism. Another framework is 
data assimilation with reduced basis, which also has extensive 
research in recent years, see [22—32]. In a word, many efforts 
in data assimilation aim to improve the accuracy, efficiency 
and robustness of physical field reconstruction. We refer the 
interested readers to the review paper [33]. 


On the other hand, the location of sensors affects the accu- 
racy and robustness of the reconstructed field, thus an impor- 
tant aspect is to study the optimization of sensor placement. 
In [22], the authors proposed the so called Generalized Em- 
pirical Interpolation Method (GEIM) [34] to select the quasi- 
optimal sensor locations in the framework of data assimila- 
tion and reduced basis. Subsequently, validation was con- 
ducted on three types of reactors in operation at Electricité de 
France (EDF). In [35], the authors applied simulated anneal- 
ing to optimize the placement of fixed in-core detectors, using 
the variance-based and information entropy-based methods to 
define the objective function. Recently, clustering, such as the 
K-means algorithm is used to optimize the in-core detector lo- 
cations for flux mapping in AHWR [36, 37]. In recent work, 
greedy algorithm is used to optimize sensor locations on a 
grid, adhering to user-defined constraints, in building the nu- 
clear digital twins based on the Transient Reactor Test facility 
(TREAT) at Idaho National Laboratory (INL) [38]. All these 
methods attempt to optimize the placement of the in-core de- 
tectors in a heuristic manner, and they are limited to a fixed 
sensor arrangement as that used in a training process; how- 
ever, little research has been done on algorithms for handling 
detector vibrations. 


The vibration of the in-core sensors near their nominal lo- 
cations is a new issue, which may be caused by the aging 
of operational reactors. A typical limitation stems from the 
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fact that all methods mentioned above fail to handle spatially 
moving sensors. Recently, the work in [39] opens a new path- 
way toward the practical use of neural networks for global 
field estimation, considering sensors could be in motion and 
could become online or offline over time. In that work, the au- 
thors used the Voronoi tessellation [40] to obtain a structured- 
grid representation from sensor locations, thereafter the con- 
volutional neural networks (CNN) can be used for building 
the map from movable sensors to the physical field. Inspired 
by that work, we adapted the framework to the application of 
field reconstruction in nuclear reactors, which is able to take 
the vibration of sensors into consideration when reconstruct- 
ing neutronic fields. 

The rest of this paper is organized as follows. In Section II, 
we give a detailed description of the methodology for field 
reconstruction with movable sensors using Voronoi tessella- 
tion in combination with convolutional neural networks (V- 
CNN). In Section III, we present the physical model and the 
detailed process to reconstruct the neutronic field. Section IV 
illustrates the numerical results, in which various error met- 
rics have been presented to evaluate the performance of the 
method. Finally, we give a brief conclusion and further works 
in Section V. 


Il. METHODOLOGY FOR FIELD RECONSTRUCTION 


WITH MOVABLE SENSORS 


Our goal is to reconstruct a two-dimensional neutronic field 
€ R?**" in the reactor core domain Q € IR? from sparse 
and limited in-core sensor observations y € IR"^^* at loca- 
tions ry,, i = 1,..., obs. Here Nops indicates the number of 
in-core sensor observations and ng and ny denote the number 
of grid points in the horizontal (x) and vertical (y) directions 
on a high-resolution field, respectively. The challenge here is 
to handle movable sensors at their nominal locations over the 
field. The field reconstruction process should be performed 
with only a single machine learning model to avoid retrain- 
ing when sensors move from their nominal locations. This is 
achieved through two key processes, i.e., 


(i) a partition method using Voronoi tessellation which is 
able to tolerate the local perturbations of sensor loca- 
tions; 


(ii) a machine learning framework that maps the observa- 
tions to the global physical field in the same structure. 


We remark here that, considering the sensors in the core 
of a reactor are fixed, such as self-power neutron detectors 
(SPND), we only consider cases of sensors vibrate near their 
fixed positions, rather than significant movement in the whole 
domain over time. The latter corresponds to the cases referred 
to [39]. 

We illustrate in Fig. 1 the framework, i.e., Voronoi 
tessellation-assisted convolutional neural networks (V-CNN) 
for neutronic fields reconstruction, and give a detailed de- 
scription of each component in the following sections. 
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A. Voronoi tessellation for spatial domain partitions 


To achieve the goal of the reconstruction using mov- 
able sensors, the Voronoi tessellation is a key step that 
maps the observations to the whole spatial domain. For a 
given space Q, which is generally in 2D, a set of points 
{ri, i = 1,..., Nobs} € Q. The tessellation approach op- 
timally partitions the given space €? into noy, regions G = 
(91,92, ..., 9n,,,] using boundaries determined by the mea- 
sure d among the given points. Using the measure d, Voronoi 
tessellation can be expressed as 

gi = ir € Q | d(r, ri) < d(r,rj), ji]. (1) 
In this article, the Euclidean measure is used and the Voronoi 
boundaries between points are their bisectors. Fig. 2 illus- 
trates an example of 81 points and the related Voronoi par- 
titions in Euclidean measure. The Voronoi tessellation pro- 
vides a convenient way to project the sparse sensor observa- 
tions to the global physical field, thus leaning the map from 
sparse observations to the global field using CNN is possible. 
More importantly, this partition process is able to tolerate the 
local perturbations of sensor locations. For more details on 
the mathematical theory of Voronoi tessellation, we refer the 
readers to the studies such as [40—43]. 


B. Input and output of the machine learning model 


To reconstruct the physical field using machine learning 
method, we utilize the following process to prepare the input 
data of model. 


(i) Determine the sensor locations r,,, ? = 1,..., Nobs in 
the reactor core. ry, may vibrate from its nominal lo- 
cation rpominal, and we have ry, = groceries 
where ór,,, is a small quantity caused by vibration. 


(ii) Calculate the Voronoi tessellation s; € IR"**""v us- 
ing ry, i = 1l,..,mog,. The Voronoi tessellation 
first partitions the reactor domain Q into no», regions 
G = (91, gnap, } With Q = U7*' gi, each region gi 


contains one sensor located at ry,, and s; is defined as 


follows 
si(r) = { 


(iii) Prepare the Voronoi mask field m € IR"**"v using 
Voronoi tessellation s;, i = 1,..., Nobs. The element 
om/(r) satisfies 


1 ifreg; 
0 otherwise 


(2) 


Qm(r) = Yr; , ifr € gi. (3) 


The Voronoi mask field ¢,,, and the related target neutronic 
field n are provided to a machine learning model F5, such 
that Fs : dm —9 $4, where ô means that the model is trained 
for a given 6. The final output of the model is then denoted 
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Voronoi field reconstruction 


Deep learning-based (CNN) field reconstruction 


Input: Voronoi mask field q® = @m(u) € R%=*%y 


High-fidelity field 
M 


oronoi mask field from sensor: 


Convolutional neural network Pr = F (bm, W) = qm» € R'*XNy 


epicrain) Predicted high-fidelity field 
J 


m 


$ (train) 


n 


Label: High-fidelity field ¢,() € RN=*%y 


Fig. 1. Voronoi tessellation-assisted convolutional neural networks for neutronic fields reconstruction from discrete sensor locations in a 
two-dimensional reactor core. The input Voronoi field is constructed from 81 sensors from each center of fuel assembly. The Voronoi field 
is then fed into a convolutional neural network with the Voronoi mask field $,,, and the output of CNN is the reconstructed field ¢,.. In the 
mask field, a grid with a sensor ? (black circle) has a value of m,i, which reflects the detected value of the underlying field at site. 


Fig. 2. An example of 81 points and the related Voronoi tessellation 
using Euclidean measure. 


by ¢, = F5(¢m), where the subscript ‘r’ presents the re- 
constructed neutronic field. With the specified input vector 
holding the observed value and the position information of 
the sensors, the proposed model can deal with arbitrary sen- 
sor locations and arbitrary numbers of sensors. 


Note here that the perturbations of sensor locations and 
whether this would have an impact on the effectiveness of 
the network is also investigated in this work. The typical 
amplitude ó of the vibration of sensors in a reactor core is 
less than 1cm [6], to go further, we investigated the cases of 
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6 = lem,3cm and 5cm, and the effects of the number and 
location of sensors in these scenarios will be illustrated in nu- 
merical results section. 

To construct the training set {¢@™, (trai) and test set 
(950. ot} for the learning process, a physical model of 
the underlying problem 


(4) 


is solved numerically. Here y € D C IR? is the p-dimensional 
parameter of the model and D is the feasible parameter do- 
main. The training set and test set are accumulated by solving 
Eq. (4) over a discrete set (discere? which is representative of 
D. 


C. Learning the map using convolutional neural networks 


Once the input data is prepared for the model, a CNN 
model can be used to learning the map from observations to 
the field, in the same way as handling images [44—46]. In this 
work, the channel of CNN is set to one, and for each layer the 
process to extract key features of input data through filtering 
operations can be expressed as 


H—1H-1 


(1-1) 
c (9 5 titp-0,j+e-0,mħpem + bm). 
p=0 c=0 


l 
ib. 


(5) 


Here C = floor(H/2), q- and q are the input and out- 
put data at layer l, respectively; hpem presents a filter of size 
H x H and bpm is the bias. lmaz, H and m denote the number 
of layers, the size of filter and the number of filters, respec- 
tively. The output of each filter operation is fed to an acti- 
vation function o(-) as the output of neurons. In this work, 
we chose the rectified linear unit (ReLU) o(z) = max(0, z) 


as the activation function [47]. We use the ADAM optimizer 
with an early stopping criterion for training, and a threefold 
cross-validation is used [48—50]. Furthermore, the Lo norm is 
used to measure the error in the learning process. The detail 
parameters of the CNN used in this work are summarized in 
Table 1. 


We fed the Voronoi mask field ¢,, € R"**"" to the model 
F, ie, qU = om, and the output of model is a high- 
resolution neutronic field qe), The learning process can 
be formulated as 


215 w = arg min ||g — F(¢m, w)||2, (6) 


where w denotes model parameters, specifically filters of 
CNN in this work. Once the training process is finished, the 
field reconstruction can be achieved by feeding the observa- 
tions y to the model F, i.e., ó(y) = F(dm(y), w). 
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Table 1. Parameter settings of CNN. 
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Fig. 3. Geometry of 2D IAEA nuclear core, upper octant: region 


Layers Hidden Filter size Number of filters Learning rate Resolution assignments, lower octant: fuel assembly identification (from [51]). 


lmaz layers H m of ADAM 


Ng X Ny 
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II. APPLICATION TO NEUTRONIC FIELD 


RECONSTRUCTIONS 


220 
221 


A. Physical model 


222 


In this section, we test the reconstruction method in nuclear 
reactor core. We consider a typical benchmark in nuclear re- 
actor physics, i.e., the 3D IAEA benchmark problem [51]. 
This benchmark was prepared by the computational bench- 
mark problems committee of the mathematics and computa- 
tion division of the American nuclear society. The reason for 
choosing this benchmark is that this benchmark is adapted 
from realistic reactors, and the geometry and the composition 
are much more complex than single region or two regions 
232 problems. Once this test is passed, the next step will be to 


233 test the method directly based on real reactor calculations 


For algorithm testing purpose, we use the 2D IAEA case 
2:5 Which represents the midplane z = 190 cm of the 3D IAEA 
23 benchmark, see [51, page 437] for a detail description. The 
2D geometry of the reactor is shown in Figure 3, where only 
one quarter is given because the rest can be inferred from 
symmetry along the x and y axes. This one quarter is denoted 
by © and it is composed of four subregions of different phys- 
ical properties: the first three subregion form the core domain 
Q4,2,3, while the fourth subregion is the reflector domain (4. 
Certain Newman boundary conditions are satisfied on the x 
and y axes considering symmetry, and the zero boundary con- 
245 dition is satisfied on the external border, see Fig. 3. 


234 


244 


The neutron fields are composed of fast flux and thermal 
247 flux, i.e., 6 = ($1, $2) are modeled by two-group neutron 
zs diffusion equation with suitable boundary conditions. The 
249 fluxes are the solution to the following eigenvalue problem 
250 (see [52, 53]). To be precise, the flux ¢ satisfies the following 
es eigenvalue problem: Find (4, 9) € C x L?*(Q) x L?*(Q), 
252 S.t. 


246 


1 
—V(DiVo1) + (Zar + E12 + DiBZ) & = xvEgaóo 


— V (D2V 92) + (Ea; + D2B?,)b2 = E1201 

253 (7) 
254 With the zero boundary condition ¢; = 0, i = 1, 2 on external 
ess border OQ and Newman boundary conditions 0(¢;)/O(n) = 
26 0, 2 = 1,2 on axises. The generated nuclear reactor rate 
27 is P = v¥iz1¢1, which reflect the power distribution over 
258 the core. The following parameters are involved in the above 
259 equation: 

* Dj: the diffusion coefficient of group i with i € {1, 2}; 


260 

261 e X44: the macroscopic absorption cross section of 

262 group 1; 

e 1-42: the macroscopic scattering cross section from 
group 1 to 2; 


263 
264 


265 e X p: the macroscopic fission cross section of group i; 


266 * v: the average number of neutrons emitted per fission. 
27 The axial buckling B2; = 0.8- 107? for all regions and energy 
28 groups. The nominal values of the coefficients in the diffusion 
ze» model (7) are listed in Section III A. 
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Table 2. Parameter values of the 2D IAEA benchmark problem. 


"E 


Region Dı D» i152 Nay Xa 2  vXyi vV%ş 2 Material 
(cm) (cm) (cm!) (em) (cm7') (cm ^?) (cm?) 
Qı 150040 0.00 0.01 0.080 0.00 0.135 Fuell $; 
Q2 150 0.40 0.02 0.01 0.085 0.00 0.135 Fuel2 -EEEEEEEE 3 
Q3 1.50 0.40 0.02 0.01 0.130 0.00 0.135 Fuel 2 + Rod: j i 
Q4 2.00 0.30 0.04 0.00 0.010 0.00 0.000 Reflector E: : | Q2 
* The axial buckling B2; = 0.8 - 1074 for all regions and energy b 


groups. 


Under some mild conditions of the parameters, the max- 
imum eigenvalue Amax is real and strictly positive (see [54, 
Chapter XXI]). The associated eigenfunction ¢ is also real 
and positive at each point x € Q and it is the flux of interest. 
In neutronics, it is customary to use the inverse of Ag, that 
is called the multiplication factor 


1 
Amax l 


Here, for each parameter setting, keff is determined by the 
solution to the eigenvalue problem (7). The maximum eigen- 
value Amax can be computed based on the well-known power 
method (see [52]). In this work, the Freefem--4 [55] is used 
to solve the 2D IAEA benchmark problem. The spatial ap- 
proximation uses P4 finite elements with a grid of size h = 1 
cm, thus the resolution of the field is n; x ny — 171 x 171. 


keff : 


(8) 


B. Field reconstruction 


To simulate the variation of the neutronic fields with re- 
spect to parameter variations, we take the parameters in Sec- 
tion III A as uncertain parameters. To be precise, we assume 
that 


HZ (m, H2; Un) eD-c [Lino minat 0.8, li nomina 1.2]" 

(9) 
where D is the parameter domain. We randomly gener- 
ate 10000 samplings of u in D, and solve Eq. (7) to get 
a collection of 10000 samples of neutronic fields M 
{¢1(), d2(u), P(x) | u € D). Among then, 8000 samples 
are used for training, 1000 samples are used for validation 
and 1000 for testing. In this work, we train three CNN mod- 
els with the same input data, i.e., the thermal flux $2 which 
is measurable with in-core detectors, and the output fields are 
the fast flux ¢,, the thermal flux $» and the reaction rate P, 
see Fig. 4 for a schematic diagram. 

To synthesize observations, we assume in the center of each 
assembly, there exists an in-core sensor to acquire the thermal 
flux. We further assume that these sensors can move in a local 
square with width 6 cm, centered at the center of each assem- 
bly. In this work, we will brought out the numerical tests for 
the cases 6 = 1,3,5 cm, to investigate the effect of differ- 
ent levels of vibration of sensors. This means that the obser- 
vations are generated randomly from the windows of width 
6 centered at their nominal locations, see Fig. 5. Then the 
model F is trained based on the set M, a schematic diagram 
for the training process can be found in Fig. 1 and Fig. 4. 


312 


324 


TTT 


Fig. 4. A schematic diagram for the reconstruction the neutronic 
fields. 


C. Error metrics 


Before we present numerical results, we define several met- 
rics to evaluate the quality of various field reconstructions. 
The normalized root-mean-square residual of the difference 
of the reconstruction @,. and the test field œ; is 


Neal 
SIE IE 


In nuclear engineering domain, the error of the reconstructed 
field in Lo is another import metric, which reflects the worst 
case for each reconstruction. 


— dés — dello 
Set e ^ 


The total average root-mean-square residual and the standard 
deviation over the given test set M are defined as 


E(ex(ó)) ^ :— averagege m (ex (9)) 
STD(e,()) := standard deviationge «(ey (¢)) ' 


(10) 


(1) 


(12) 


where x denotes La or Loo norm. 
Furthermore, the average assembly field (fluxes and power 
rate) and the related errors are also investigated. The average 


assembly power is defined as 
1 
] eae, 
Uk Ju, 


where ¢ denotes $1, ¢2 or P, vj, denotes the volume of the 
k-th subassembly, and k designates the fuel assemblies as 
shown in lower octant of Fig. 3. The average error e2(@ass) 
and maximum relative error €..(@ass) and the related aver- 
age E(e,(@ass)) and standard deviation ST D(e,(¢ass)) of 
the errors over the given test set M can also be defined in a 
similar way. 

To this end, we introduce the structural similarity (SSIM) 
[56] index to measure the field reconstruction. Contrary to 
the general Lo error, the SSIM index provides a measure of 
the similarity by comparing two images based on luminance 
similarity, contrast similarity and structural similarity infor- 
mation. 


Óass,k = (13) 
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Fig. 5. Different amplitudes of vibrations of a sensor near the nominal location in an assembly. 


IV. NUMERICAL RESULTS 


We have finished to describe the methodology for field re- 
construction with movable sensors using CNN in Section II 
and presented a detailed process for neutronic fields recon- 
struction based on a typical benchmark nuclear engineering 
domain in Section III. In this section, we illustrate the numer- 
ical performance of the proposed method. 


A. Performance for the benchmark problem 


We first show in Fig. 6 the error distribution of the recon- 
structed fields for different vibration amplitudes, i.e., ô = 
1,3, 5 cm for the 2D IAEA benchmark problem. The recon- 
struction of $4 using observations from $» shows the best 
performance than the reconstruction of $9 and P. Further- 
more, with the amplitude of the vibration becoming larger, 
i.e., 6 varies from 1 cm to 5 cm, the reconstructed error also 
increases. The largest error appears around the interface of 
the fuel and reflector, because of the discontinues of materi- 
als, particularly for the fields of thermal flux @2 and the power 
rate P. 

The averaged assembly values of the reconstructed fluxes 
and power rate are shown in Fig. 7, with which the same con- 
clusions can be drawn. Note that because of the average pro- 
cess, the relative errors in assembly wise are much smaller 
than that of the pin wise. 

More importantly, three main conclusions can be drawn 
from the prior analysis of the numerical results: 


(i) The proposed V-CNN is able to reconstruct the multi- 
field with observations only from thermal flux; 


(ii) The reconstruction errors in assembly wise are far be- 
low 596, which is acceptable for engineering usage, i.e., 
less than 10%, which is a normal criteria in reactor 
physics, see [57] for more information; 


(iii) Even with a movement of amplitude 6 = 5 cm for 
sensors, the proposed V-CNN is able to reconstruct the 
field with an error less than 1046. 
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B. Average performance over a test set 


In order to investigate the generalization ability of the field 
reconstruction method, we analyzed the error performance 
on 1000 test samples. The error metrics are average relative 
Lo error E'(e2(¢)), average relative La error E'(e..(¢)), av- 
erage relative assembly Lo error E(e2(¢ass)), average rela- 
tive assembly La error E (eso(Pass)) and the average SSIM, 
E(SSIM(Q)), and the standard deviation of the afore men- 
tioned errors. 


Table 3 illustrates the numerical results of the errors for the 
reconstruction of $, over the 1000 test samples. All the error 
metrics present a good agreement between reconstructed field 
and the original field over the test set. Notice that the worst 
errors, i.e., the Loo errors both in pin wise and assembly wise 
are below 2%. This good performance is attributed to the 
smoothness of the fast flux. This is because the fast flux has 
a relatively longer diffusion length than that of thermal flux, 
thus the fast flux is less affected by the heterogeneous of ma- 
terials of this benchmark, see Fig. 6(a) for example. 


Table 4 and Table 5 illustrate the numerical results of the 
errors for the reconstruction of @2 and P over the 1000 test 
samples. The average Lo errors in pin wise of the thermal 
flux and the power rate are below 10% when the vibration 
amplitude is less than 3 cm. Once the amplitude of the vibra- 
tion gets larger, we will get an average Loo error larger than 
10%, which in some sense is not acceptable in practical engi- 
neering application. On the other hand, if we look the errors 
in assembly wise, the Loo errors are much smaller. The worst 
case appears for 6 = 5 cm when reconstruct $5, which leads 
to an error of E(es;($455)) = 0.0288 with standard deviation 
ST D(es(Pass)) = 0.0109. Again, this results confirm the 
acceptance for engineering application. Notice that though 
the relative Lo errors amount to 10%, most of the points ap- 
pear around the interface between the fuel and the reflector. 
The relative large error in this domain is not crucial for safety 
analysis. To this end, the SSIM indexes in all the three ta- 
bles are larger than 0.99, which again demonstrate excellent 
performance for all the field reconstructions. 
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Fig. 6. The reconstructed fields for different vibration amplitudes, ô = 1, 3, 5 cm for the 2D IAEA benchmark problem. 
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Fig. 7. The reconstructed fields in assembly wise for different vibration amplitudes, 6 = 1,3,5 cm for the 2D IAEA benchmark problem. 
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Table 3. The errors in different metrics for the reconstruction of $1 
using thermal flux observations from movable sensors. 


Width 1 3 5 


E(e2(¢)) 0.0097 0.0103 0.0119 
ST D(e2($)) 0.0017 0.0020 0.0031 
E(ess(d)) 0.0276 0.0293 0.0303 
ST D(ess($)) 0.0078 0.0090 0.0091 
E(e€2(¢ass)) 0.0084 0.0094 0.0105 
ST D(e2(bass)) 0.0020 0.0022 0.0032 
E(esc(Qass)) 0.0155 0.0159 0.0171 
ST D(es(Pass)) 0.0046 0.0054 0.0064 
E(SSIM(Q)) 0.9986 0.9982 0.9979 
ST D(SSIM(¢)) 0.0003 0.0004 0.0004 


Table 4. The errors in different metrics for the reconstruction of $2 
using thermal flux observations from movable sensors. 


Width 1 3 5 


E(e2(¢)) 0.0167 0.0190 0.0213 
ST D(e2(¢)) 0.0025 0.0032 0.0041 
E(€c0(#)) 0.0751 0.1001 0.1219 
STD(eæ(¢)) 0.0197 0.0293 0.0391 
E(€2(dass)) 0.0120 0.0134 0.0141 
ST D(e2(¢ass)) 0.0025 0.0027 0.0028 
E(€co(dass)) 0.0216 0.0256 0.0288 
ST D(€.0(Pass)) 0.0085 0.0093 0.0109 
E(SSIM(¢)) 0.9969 0.9964 0.9957 
STD(SSIM(¢)) 0.0006 0.0008 0.0010 


C. Robustness analysis 


The robustness of the reconstruction with respect to the 
number of observations ns, and the amount of training data 
Nsnapshot 15 examined. We show in Fig. 8 the dependence of 
the relative reconstruction errors in Lo norm and Leo norm 
ON Nsnapshot. = 128, 1280, 4096, 8192, 15743 and on Nobs = 
25, 45,81 for recovering thermal flux $2 over the test set. 
These numbers of observations nop, = 25, 45, 81 correspond 
to a sparsity of 0.085596, 0.154%, 0.277% against the num- 
ber of grid points over the field. It can be seen from the two 
figures that the proposed method shows great robustness with 
respect to the sparsity of sensors and training data. Both low 
number of observations and low number of training data lead 
to low-level reconstruction error. Furthermore, the addition 
of training data improves the reconstruction accuracy much 
more than the addition of sensors. This result shows that the 
proposed field reconstruction framework is tolerant of sensor 
failures, and confirms the potential of real engineering appli- 
cation. 

To investigate the robustness of the recovery with respect 
to the observation noise, we added a noise €, randomly sam- 
pled in the range (—6,6) to each clean observation y, thus 
we have the noisy observation y^ = y(1 + €s) for each 
sensor. The dependence of the relative reconstruction er- 
rors in Ly norm and Leo norm of different noise level, i.e., 
o = 0.01, 0.02, 0.03, 0.04, 0.05 for the recovering of thermal 
flux #2 are shown in Table 6. The test is carried out with 81 
sensors with vibration amplitude 6 = 5 cm. The errors are 


Table 5. The errors in different metrics for the reconstruction of P 
using thermal flux observations from movable sensors. 


Width 1 3 5 


E(e2(@)) 0.0137 0.0164 0.0253 
ST D(e2(¢)) 0.0031 0.0046 0.0060 
E(e.0(¢)) 0.1429 0.2069 0.2640 
STD(es(ó) 0.0811 0.1079 0.1119 
E(€2(¢ass)) 0.0097 0.0108 0.0192 
ST D(e2(¢ass)) 0.0027 0.0037 0.0058 
E(€co(¢ass)) 0.0182 0.0191 0.0257 
STD(€cc(bass)) 0.0069 0.0081 0.0082 
E(SSIM(Q)) 0.9951 0.9947 0.9920 


STD(SSIM(@)) 0.0012 0.0018 0.0023 


444 firstly averaged over 100 repeated random observation sam- 
44s plings for each field reconstruction, and then averaged over 
4s the test set. On average, there is a significant change in the 
447 reconstruction error once the observation is polluted by noise. 
4 The reconstruction error shows a slow linear growth trend 
4 With respect to the noise level. Though the Lə error is be- 
450 low 10%, which is satisfactory for nuclear engineering appli- 
cations, the Lo error still stands around 30%, which is not 
so satisfactory, provides a direction for further research in di- 
minishing the Loo error. 


45 


452 


453 


Table 6. The dependence of the relative reconstruction errors in L2 
norm and Loo norm of noise level c for the recovering of thermal 
flux $2. 


Noise level E(e2($)) E(ess(Q)) 


0.0 0.0213 0.1219 
0.01 0.0493 0.2956 
0.02 0.0510 0.2975 
0.03 0.0536 0.3013 
0.04 0.0572 0.3078 
0.05 0.0616 0.3158 


V. CONCLUSIONS 


454 


45 X In this article, a Voronoi tessellation-assisted convolutional 
neural networks (V-CNN) is proposed for neutronic fields re- 
construction to settle the vibrations of in-core sensors, which 
may be caused by the aging of operational reactors. The ob- 
servations from movable in-core sensors are projected to the 
same global field structure, this projection is achieved with 
Voronoi tessellation, holding the magnitude and location in- 
formation of sensors. The general convolutional neural net- 
455 Works were used to learn the map from observations to the 
4&4 global field. Furthermore, the proposed method is able to re- 
4s construct the multi-physics fields e.g., the fast flux, thermal 
ass flux and power rate distributions using observations from sin- 
47 gle field e.g., thermal flux. 

Numerical tests based on IAEA benchmark proved its effi- 
s9 ciency of the proposed method. Three main conclusions can 
470 be drawn from the prior analysis of the numerical results: 
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Fig. 8. The dependence of the relative reconstruction errors in L2 norm and Loo norm on the number of training snapshots Nsnapshot = 
128, 1280, 4096, 8192, 15743 and the number of observations nos, = 25, 45, 81 for recovering thermal flux ¢2. 


servations only from thermal flux; 


(ii) All the reconstruction errors in average are below 5%, 
which is totally acceptable for engineering usage; 


(iii) Even with a vibration amplitude of 6 = 5 cm for sen- 
sors, V-CNN still performs well. 


The original CNN framework is used for image processing, 
which is now adapted for field reconstruction with rectangu- 
lar mesh division in this work. For field reconstruction with 
irregular mesh, an additional mesh mapping is necessary to 
map the irregular mesh to a rectangular mesh. The adaptabil- 
ity of the proposed method to various reactor configurations 
would be a continuous work of this article. 

This prior study provides a new approach for dealing with 
field reconstructions with vibration sensors. Future works 
are possibly to bring out the uncertainty quantification of 
V-CNN considering observation noise systematically and to 
push forward to the practical engineering applications based 
on real nuclear reactors e.g., the HPR1000 reactor developed 
in China [58]. In this aspect, to evaluate the data uncertainty, 
the probabilistic neural network [59] or Bayesian neural net- 
work [60] could be investigated for a combination of V-CNN; 
to evaluate the epistemic uncertainty of the model, the Gaus- 
sian stochastic weight averaging technique [61] or other tech- 
niques could be investigated. 

To investigate the adaptability of the proposed method to 
the HPR1000 reactor, a pin-by-pin wise field calculation is 
necessary to consider the fuel and sensor vibration, which is 
now in the process of our group. However, in practical en- 
gineering case, the vibrations of reactor components such as 
fuel and in-core sensors lead to very complex phenomena in 
the core. Many works [2-6] have been brought to analyze the 
induced variation of the neutronic fields (also called neutron 
noise), considering the induced variation of cross-section pa- 
rameters of the neutron diffusion equations. Inspired by the 
process of neutron noise analysis, a synthetic modeling ap- 
proach is also necessary to consider the effects of component 
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vibration. This approach is useful to clarify the interplay or 
distinctions between the field reconstruction with in-core sen- 
sor vibrations and the general reactor noise analysis. 

In addition, the combination of V-CNN and fault diagno- 
sis [62] is also a possible research point in the future. With 
the development of machine learning in the field of nuclear 
physics [63, 64], adaptions of V-CNN to nuclear physics 
where CNNs are used [65—67] are also worth trying. 
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